function [F, G] = exp_shares_NH_CES(outputs, inputs)

p_a = inputs(1);
p_m = inputs(2);
p_s = inputs(3);
EXP = inputs(4);
epsilon_a = inputs(5);
epsilon_m = inputs(6);
epsilon_s = inputs(7);
sigma = inputs(8);
omega_a = inputs(9);
omega_m = inputs(10);
omega_s = inputs(11);

% outputs(1) = P
% outputs(2) = C
% outputs(3) = C_a
% outputs(4) = C_m 
% outputs(5) = C_s
% outputs(6) = omega_a
% outputs(7) = omega_m
% outputs(8) = omega_s

G(1) = omega_a * (p_a/outputs(1))^(1-sigma) * (EXP/outputs(1))^(epsilon_a-1+sigma) - (outputs(3)*p_a)/EXP;
G(2) = omega_m * (p_m/outputs(1))^(1-sigma) * (EXP/outputs(1))^(epsilon_m-1+sigma) - (outputs(4)*p_m)/EXP;
G(3) = omega_s * (p_s/outputs(1))^(1-sigma) * (EXP/outputs(1))^(epsilon_s-1+sigma) - (outputs(5)*p_s)/EXP;
G(4) = p_a*outputs(3) + p_m*outputs(4) + p_s*outputs(5) - EXP;
G(5) = (omega_a * outputs(2) ^ (epsilon_a) * p_a ^ (1-sigma) + omega_m * outputs(2) ^ (epsilon_m) * p_m ^ (1-sigma) + ...
    omega_s * outputs(2) ^ (epsilon_s) * p_s ^ (1-sigma))^(1/(1-sigma)) - EXP;

%G = [G1, G2, G3, G4, G5];
F = [];

